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The correct sizing of the heat transfer surfaces is important to ensure proper operation, load turndown, 
and optimization of circulating fluidized beds (CFBs). From this point of view, in this study, the thermo¬ 
dynamic second law analysis of heat transfer surfaces in CFBs is investigated theoretically in order to 
define the parameters that affect the system efficiency. Using a previously developed 2D CFB model 
which uses the particle-based approach and integrates and simultaneously predicts the hydrodynamics 
and combustion aspects, second law efficiency and entropy generation values are obtained at different 
height and volume ratios of the heat transfer surfaces for CFBs. Besides that, the influences of the water 
flow rates and heat exchanger tube diameters on the second law efficiency are investigated. Through this 
analysis, the dimensions, arrangement and type of the heat transfer surfaces which achieve maximum 
efficiency are obtained. 

© 2008 Elsevier Ltd. All rights reserved. 


1. Introduction 

Superior environmental performance of CFBs is one of the prime 
motivations of its extensive use in industry. CFBs are also charac¬ 
terized by their approximate isothermal nature and high rates of 
heat transfer between the fluidized medium and the heat transfer 
surfaces. Design or scale-up of heat transfer systems in CFBs re¬ 
quire estimation of the effective heat transfer coefficient on heat 
transfer surfaces in contact with the fluidized medium. It is there¬ 
fore essential to understand the mechanisms of heat transfer in 
CFBs, and to develop an appropriate model to predict the rate of 
heat transfer. 

There have been many studies concerning the fundamental 
analysis of the heat transfer between the fluidized beds and heat 
transfer surfaces in both laboratory and industrial-scale units. Basu 
and Nag [1] and Glicksman [2] presented comprehensive reviews 
of CFB heat transfer. Recently, Grace et al. [2] experimentally and 
theoretically investigated the effects of various operating parame¬ 
ters on the heat transfer process. Xie et al. [3] proposed a 2D model 
which coupled radiation, conduction and convection from the hot 
core on the furnace side to conduction and convection into the 
coolant on the wall side. In a further study authors also developed 
an advanced model to accommodate the 3D membrane wall geom¬ 
etry Xie et al. [3]. Gupta and Reddy [4] proposed a mechanistic 
model to predict the bed-to-wall heat transfer coefficient in the 
top region of a CFB riser column for different riser exit configura¬ 
tions. Chen et al. [5] suggested that mechanistic models based on 
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the surface renewal concept hold promise for design and scale- 
up of heat transfer systems for both bubbling dense beds and fast 
CFBs. 

Although there have been many studies concerning the funda¬ 
mental analysis of the heat transfer between the fluidized beds 
and their heat transfer surfaces as mentioned above, the second 
law analysis of the CFBs is very rare in the literature [6-12]. It is 
a well known fact that thermodynamics is concerned not only with 
the conservation of energy but also with quality of energy. The tra¬ 
ditional first law analysis based upon component performance 
characteristics together with energy balance can always lead to 
the correct final result. But, this analysis cannot locate and mea¬ 
sure the source of losses. This is mainly because the first law 
embodies no distinction between work and heat, no provision for 
quantifying the heat and no accounting for work lost. The first 
law efficiency does not take into account the quality of energy. It 
has long been understood that traditional first law analysis is 
needed for modeling the engine processes but a novel approach 
is also required to give the engineer the best insight into the en¬ 
gine’s operation. The heat transfer surfaces directly affect the 
hydrodynamic behavior and thus fluidized beds’ combustion per¬ 
formance. As for studies concerning fluidized beds’ heat transfer 
surfaces, there are only two investigations taking into account 
the exergy analysis of heat transfer surfaces. The first study con¬ 
ducted by Eskin and Kilic [8], investigates the estimation of cooling 
tube location in bubbling fluidized bed coal combustors through 
exergy analysis. Whereas the second study conducted by Gungor 
and Eskin [10] takes into account the thermodynamic analysis of 
heat transfer to the immersed surfaces in small-scale CFBs. It is 
clearly observed that the effect of heat transfer surfaces on the 
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Nomenclature 



A 

area, m 2 

t 

time, s 

Ar 

Archimedes number 

u 

overall heat transfer coefficient, kW/m 2 I< 

a 

availability, kj/kg 

Uo 

superficial velocity, m/s 

C 

gas concentration, kmol/m 3 

Ut 

particle terminal velocity, m/s 

c 

specific heat capacity, kj/kgK 

X k 

char mass fraction, kg-char/kg-bed material 

Dg 

dhe 

diffusivity coefficient for oxygen in nitrogen, m 2 /s 
heat exchanger tube diameter, m 

z 

height above the distributor, m 

d p 

particle diameter, m 

Greek symbols 

dpi 

particle dimension interval, m 

X 

available energy, exergy, kW 

e 

emissivity 

A T m 

log-mean temperature difference, °C 

g 

gravity, m/s 2 

S 

thickness of the annulus, m 

H 

specific enthalpy, kj/kg 

£ 

void fraction 

H b 

bed height, m 

£ c 

averaged fraction of wall area covered by clusters 

Ff U pp 

upper zone height, m 

Op 

volumetric concentration of particles in the dispersed 

h 

overall bed-to-wall heat transfer coefficient in the upper 


phase 


zone, W/m 2 K 

0 

height ratio 

h B 

overall bed-to-wall heat transfer coefficient in the bot¬ 

y 

volume ratio 


tom zone, W/m 2 K 

mi 

second law efficiency 

^HE 

starting height of the heat transfer surface (m) 

p 

density, kg/m 3 

h g 

convection heat transfer coefficient for gas phase, W/m 2 

K 

convection heat transfer coefficient for particle phase, 

o 

Stefan-Boltzman constant, W/m 2 I< 4 

h p 

Subscripts 


W/m 2 K 

b 

bed 

hr 

radiative heat transfer coefficient, W/m 2 K 

g 

gas 

k c 

char combustion reaction rate, kg/s 

he 

heat transfer surfaces 

k g 

gas conduction heat transfer coefficient, W/m K 

i 

Compartment number 

m 

mass, kg 

in 

input 

m 

mass flow rate, kg/s 

L 

lost 

Pr 

Prandtl number 

0 

reference 

a 

amount of heat transferred, kW 

out 

output 

R 

combustor radius, m 

P 

particle 

Rhe 

the ratio of overall heat exchanger area to heat exchan¬ 

t 

total 


ger length, m 2 /m 

w 

water 

SH 

Sgen 

T 

Sherwood number 
entropy generation, kW/K 
temperature, K 

wall 

wall 


CFB combustion efficiency remains as a domain to be investigated 
in detail. This manuscript deals with the effect of the heat transfer 
surfaces on the dimensions, arrangement and types of surfaces in 
order to achieve maximum efficiency which is the main contribu¬ 
tion of the manuscript to the literature. This study proves that both 
the location and the dimensions of the heat exchangers are extre¬ 
mely important design parameters of the CFB combustors for opti¬ 
mization of system efficiency. 

As for the thermodynamics, exergy is defined as the maximum 
amount of work which can be produced by a system or a flow of 
matter or energy as it comes to equilibrium with a reference envi¬ 
ronment. Unlike energy, exergy is not subject to a conservation law 
(except for ideal, or reversible, processes). Rather exergy is con¬ 
sumed or destroyed, due to irreversibilities in any real process. 
Second law (exergy analysis) quantifies and locates the losses that 
help in optimizing the thermal systems. In order to evaluate the 
inefficiencies associated with the various processes - second law 
analysis must be applied [13]. For second law analysis, the key con¬ 
cept is “availability”’ (or exergy). The availability of a material is its 
potential to do useful work. Unlike energy, availability can be de¬ 
stroyed with such phenomena as combustion, friction or mixing. 
The reduction of irreversibilities can lead to better engine perfor¬ 
mance through a more efficient exploitation of fuel, better environ¬ 
mental impact in general and higher potential savings. The use of 
exergy principles improves the engineers’ understanding of ther¬ 
mal and chemical processes and allows sources of inefficiency to 


be quantified. From this point of view, in this study, the thermody¬ 
namic second law analysis of heat transfer surfaces in CFBs is ana¬ 
lyzed in order to define the parameters that affect the system 
efficiency. Using a previously developed 2D CFB model [14,15] 
which uses the particle-based approach and integrates and simul¬ 
taneously predicts the hydrodynamics and combustion aspects, 
second law efficiency and entropy generation values are obtained 
at different height and volume ratios of the heat transfer surfaces 
for CFBs. Besides that, the influences of the water flow rates and 
heat exchanger tube diameters on the second law efficiency are 
investigated. It must be noted that since the nonexistence of exper¬ 
imental data in the literature, in this study only theoretical analysis 
is given. 

2. System analysis 

The use of CFB modeling enables the analysis of a combustion 
system involving fluid flow, heat transfer, and combustion and pol¬ 
lutant emissions. The thermodynamic analysis of heat transfer in 
circulating fluidized beds is very important in defining the param¬ 
eters that affect the system efficiency. For this purpose, thermody¬ 
namic analysis of the heat exchanger surfaces in the CFB 
combustor is evaluated through a model, which can be employed 
to simulate a wide range of operating conditions [14,15]. The 
scheme of the CFB used in the model is shown in Fig. 1. The model 
addressed in this paper uses particle-based approach which con- 
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Flue Gases 



Fig. 1 . The scheme of the CFB. 


siders two-dimensional motion of single particles through fluids. 
According to the axial solid volume concentration profile, the riser 
is axially divided into the bottom zone and the upper zone. The 
bottom zone in turbulent fluidization regime is modeled in detail 
as two-phase flow which is subdivided into a solid-free bubble 
phase and a solid-laden emulsion phase. A single-phase back-flow 
cell model is used to represent the solid mixing in the bottom zone. 
In the upper zone core-annulus solids flow structure is established. 
It is assumed that the particles move upward in axially and move 
from core to the annulus region in radially. In the annulus region, 
the particle has a zero normal velocity. 

Hydrodynamic model takes into account the axial and radial 
distribution of voidage and velocity, for gas and solid phase, pres¬ 
sure drop for gas phase and solids volume fraction and particle size 
distribution for solid phase. In order to determine the validity of 


the hydrodynamic model, the model results are compared with 
test results using the same input variables in the tests as the sim¬ 
ulation program input. Hydrodynamic model results are compared 
with experimental results obtained from various CFB test rigs at 
different size in the literature which ranges are as follows: bed 
diameter from 0.05 to 0.418 m, bed height from 5 to 18 m, mean 
particle diameter from 67 to 520 pm, particle density from 1398 
to 2620 kg/m 3 , mass fluxes from 21.3 to 300 kg/m 2 s and gas super¬ 
ficial velocities from 2.52 to 9.1 m/s. The model compares with 
experimental results for void fraction, solids volume fraction, and 
particle velocity along the bed height and bed radius and different 
bed operational parameters prove the model hydrodynamic struc¬ 
ture validation axially and radially. Additionally inward solids 
mass flux along the bed radius and pressure gradient along the 
bed height is validated [16]. 



















































A. Gungor/Applied Energy 86 (2009) 1344-1353 


1347 


The combustor model takes into account the devolatilization of 
coal, and subsequent combustion of volatiles followed by residual 
char. In the model, volatiles are entering the combustor with the 
fed coal particles. It is assumed that the volatiles are released in 
emulsion phase in the bottom zone of the CFBC at a rate propor¬ 
tional to the solid mixing rate. The degree of devolatilization and 
its rate increase with increasing temperature [17]. The composi¬ 
tion of the products of devolatilization in weight fractions is esti¬ 
mated from the correlations proposed by Loison and Chauvin [18]. 

The bed material in the combustor consists of coal, inert parti¬ 
cles and limestone. The properties and size distribution of particles 
have significant influence on the hydrodynamics and combustion 
behavior in the CFBC [19]. The model also considers the particle 
size distribution due to fuel particle fragmentation [20], char com¬ 
bustion [21] and particle attrition [22]. Particles in the model are 
divided into ten size groups in the model. The Sauter mean diam¬ 
eter is adopted as average particle size. Particles in the bottom 
zone include particles coming from the solid feed and recirculated 
particles from the separator. 

Kinetics of char combustion is modeled with a shrinking core 
with attiring shell i.e., the dual shrinking core model (assuming 
that the ash separated once formed) with mixed control by chem¬ 
ical reaction and gas film diffusion. The rate at which particles of 
size Vi shrink as follows [21]; 

r(r.) = = _ 1 ^ 0)2 _ /i\ 

Kl) dr pX k) j(l//< C) j + dpj/5/i/Dg) 

The term C 02 indicates the effective oxygen concentration seen by 
the char particles burning at any point of the combustion chamber. 
The kinetic constants for the fuel used in the model are determined 
by [23,24]. The fuel properties are given in Table 1. 

It is crucial to well evaluate the mechanism of NO x formation to 
reduce NO x in the combustor. NO x formations in combustion pro¬ 
cesses result from a combination of a thermal generation process 
and fuel nitrogen oxidation. The chemical reactions with their cor¬ 
responding reaction rates for NO x emissions formation and reten¬ 
tion in the model are available in the literature [14,15]. 

CFBs offer the possibility of removing the sulfur dioxide during 
combustion by adding limestone (CaC0 3 ) directly to the combus¬ 
tion chamber. Using such limestone directly with a fuel inserted 
to the combustion chamber leads to highly efficient desulphuriza¬ 
tion. The chemical reaction with their corresponding reaction rate 
for S0 2 retention regarding the gas temperature and particle size is 
available in literature [14,15]. 

Coal combustion model calculates the axial and radial distribu¬ 
tion of voidage, velocity, particle size distribution, pressure drop, 
gas emissions and temperature at each time interval for gas and so¬ 
lid phase both for dense bed and for riser. The model has been val¬ 
idated against two different size CFBCs, which use different kinds 
of Turkish lignite, the pilot-scale 50 kW CFBC using Tuncbilek lig¬ 
nite and an industrial-scale 160 MW CFBC using Can lignite. The 


technical parameters of the industrial-scale CFBC are steam capac¬ 
ity of 485 t/h, superheated steam temperature and pressure of 
543 °C and 17.5 MPa, respectively [14]. 

Regarding the heat transfer approach, mathematical modeling 
is an attractive solution which allows developing an appropriate 
model to predict bed-to-wall heat transfer characteristics for bet¬ 
ter understanding heat transfer mechanism to enhance combus¬ 
tion performance in a much shorter time period and at lower 
costs. The gas and solid concentrations at the riser wall along the 
bottom zone can be estimated with sufficient accuracy using the 
developed model and as a function of the local bed density, the 
bed-to-wall heat transfer coefficient is calculated as the expression 
given by Basu and Nag [1] which was validated previously [25]; 

h b =40(p b )' /2 (2) 

where p b is the local bed density; p h = p( 1 - e) + Cs. In this equa¬ 
tion, C is the gas concentration and is expressed as kg/m 3 . The 
bed-to-wall heat transfer coefficient in the upper zone of CFB con¬ 
sists of the contributions of four components i.e., conduction and 
radiation components of particle phase and convective and radia¬ 
tion components of gas phase. Any part of the heat transfer surface 
comes in contact with the solid phase and the gas phase with re¬ 
spect to the rates of their voidage at the region in vicinity of the heat 
transfer surface. This assumption is based on the ability of devel¬ 
oped hydrodynamic model to make the calculations for any given 
time step [16]. The local solids volume fraction (e p ) and bed voidage 
(e) at the riser wall along the upper zone can be estimated with suf¬ 
ficient accuracy using the developed hydrodynamic model, the bed- 
to-wall heat transfer coefficient, ft, is calculated as the sum of the 
particle conductive, gas convective and particle and gas radiative 
heat transfer coefficients in the model as follows; 

h = £ p ■ h p + e ■ h g + £ p • h T , P + e • h r , s (3) 

where h p is the conductive heat transfer coefficient for solids phase 
and h g is the convective heat transfer coefficient for the gas phase. 
Heat transfer by particle conduction h p refers to the energy transfer 
due to continuous particle motion between heat transfer surface 
and inner region of the CFB and it is calculated in the model as fol¬ 
lows [26]; 

hp = ^ • 0.009 x p r 033 Ar 0 5 (4) 

a p 

Heat transfer by gas convection h g becomes significant in fluidized 
beds at large Archimedes numbers or low solids concentrations 
[27]. The convection heat transfer from the gas phase to the wall 
is estimated by the modified equation of Wen and Miller [28], 
which was given by Basu et al. [29] as 



Table 1 

Combustor dimensions and operating conditions 


Operating parameters 

Values 

Bed diameter 

0.20 m 

Bed height 

6.5 m 

Secondary air injection height from distributor 

1.5 m 

Inlet pressure (atm) 

1.12 

Coal feed rate (range) 

10-12 kg/h 

Operation velocity (range) 

2.07-5 m/s 

Bed temperature (range) 

850-900 °C 

Primary/secondary air ratio 

2/3 

Excess air (%) 

20-30 

Size of coal feed (range) 

0.0125-6.5 mm 

Mean size of sorbent feed 

0.1-0.15 mm 

Heat exchanger tube diameter (range) 

1.905-3.81 cm 

Cooling water mass flow rate (range) 

0.88-8.8 kg/s 


The radiation heat exchange is calculated similar to parallel planes 
approach in radiation [1]. The particle radiation component of the 
heat transfer coefficient is estimated using the following equation: 


r,p {( e -i- e -i)--i}(r p -r w ) 


( 6 ) 


where e p is the emissivity of the particle and e w is the emissivity of 
the heat transfer surfaces which is the property of the surface [29]. 
The gas phase radiation heat transfer coefficient from bed-to-wall is 
estimated using the following equation [29]: 




^(r 4 -rt) 

{(e-'-e-’)-l}(T-T w ) 


( 7 ) 
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Fig. 2. The volume ratio, y. 


Basu and Nag [1] reported that the emissivity of cluster or the dis¬ 
persed phase is between 0.85 and 0.95, and Wirth [30] also used the 
emissivity of 0.91 for the analysis of commercial CFB coal combus¬ 
tors. The particles are assumed to constitute a continuous absorb¬ 
ing, emitting and scattering medium in the model. The wall 
emissivity of 0.85 and the particle emissivity of 0.90 are used in 
the simulations. 

In the model, the last control volume in the annulus region 
which is in vicinity of the wall is considered as the basis of heat 
transfer coefficient calculations. The mass flow rate throughout 
this control volume can be expressed in terms of substances enter¬ 
ing and leaving the control volume, m in and rh out , as follows: 


tHin.i-l — tHout ,i — 


drrii 
d t 


( 8 ) 


The energy balance equation in this cell can be expressed in terms 
of rate of change of energy as 


thin,i-lf^in,i-l thout,i'f^out,i T Q.release,i Q.w,i Q.wall,i — 


d t 


0 ) 


where (Release,* is the total of heat derived from the combustion of 
the coal and the volatiles in the cell, Q w l , the amount of heat trans¬ 
ferred to the cooling water, is 


Q.w,i — Uh e iAh e iAT m i 

where AT m is log-mean temperature difference. 


( 10 ) 


The lost available work is directly proportional to the entropy 
generation in a non-equilibrium phenomenon of exchange of en¬ 
ergy and momentum within the fluid and at the solid boundaries. 
For the purpose of evaluating the degree of irreversibility of the 
combustor, the entropy generation expression for each control vol¬ 
ume is obtained [10]. 

^gen, t,i — Sgen(chem. react.),i T *5gen(heat trans.),i “F ^gen(fluid fric.),i (H) 

Here the first term on the right-hand side shows the entropy gener¬ 
ated due to chemical reaction, while the entropy generated due to 
heat transfer, is shown by the second term. The third term is the en¬ 
tropy generation due to fluid friction in the combustor. 

Since the model provides a way to calculate the local velocity, 
temperature and concentration distributions within the compart¬ 
ment for each phase, a detailed calculation of the available energy 
loss, for various cooling tube arrangements, is possible. Therefore, 
for each compartment, the second law efficiency is calculated in 
terms of available energies as 

riu = 1 -^ 02 ) 

Ain 

Here the lost available energy, which is the result of the irreversibil¬ 
ities in the control volume, is expressed in terms of substances 
entering and leaving the control volume m n and m k as follows; 

Xl = J2 _ J2 + Qi f i - y) 0 3 ) 

n k \ 1 / 


Heat transfer 



Fig. 3. Particle flux structure, control volumes in the upper zone and computational conditions in the CFB riser. (See above-mentioned reference for further information). 
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Fig. 4. Flow chart of the numerical solution of the model. 


where Q, is the total heat transferred in the control volume. 

The model allows dividing the calculation domain into m x n 
control volumes, in the radial and the axial directions and in the 
core and the annulus regions, respectively. The simulation geome¬ 
try of the control volumes is shown in Fig. 3. In this study the cal¬ 
culation domain is divided into 10 x 100 control volumes in the 
radial and the axial directions and in the core and the annulus re¬ 
gions, respectively. With the cylindrical system of coordinates, a 
symmetry boundary condition is assumed at the column axis. At 
the walls, a partial slip condition is assumed for the solid and the 
gas phases [14,15]. All particles (fuel and inert material) are as¬ 
sumed to be spherical. Flow chart of the numerical solution of 
the model is shown in Fig. 4. 

Inputs for the model are combustor dimensions and construc¬ 
tion specifications (insulation thickness and materials, etc.), pri¬ 
mary and secondary air flow rates; fuel feed rate and particle 
size distribution, fuel properties, inlet air pressure and tempera¬ 
ture, ambient temperature and the superficial velocity. A continu¬ 
ity condition is used for the gas phase at the top of the cyclone. The 
cyclone is considered to have 98% collection efficiency. 

3. Simulation results and discussion 

Using the presented heat transfer model, the second law effi¬ 
ciency of the CFB is obtained at various conditions. In performing 
the simulations, the volume ratio, y, which is the ratio of the vol¬ 


ume of heat exchanger in the control volume to the total heat ex¬ 
changer volume (Fig. 2); and the height ratio, 0, which is the ratio 
of the starting height of the heat transfer surfaces to the total 
height of the CFB are defined (Fig. 1). Simulation results are ob¬ 
tained from 300 kW CFB unit of 20 cm, i.d., 650 cm tall main col¬ 
umn (riser). Combustor dimensions and operating conditions are 
given in Table 1. More details of the CFB unit are given in the liter¬ 
ature [32]. The compositions of coal used in the simulations are gi¬ 
ven in Table 2. 

The influence of the positioning of the cooling tubes ( 0 ) on 
combustor efficiency is shown in Fig. 5. It plots the variations of 
combustor second law efficiency with dimensionless combustor 
height at constant volume ratio, y = 0.30 for different height ratios 


Table 2 

Proximate and ultimate analysis of fuel (Tuncbilek lignite) 


Proximate analysis (wt%) 

Moisture 

20.80 

Ash 

23.70 

Volatile matter 

27.50 

Fixed carbon 

41.30 

Ultimate analysis (wt%, dry) 

C 

59.29 

H 

4.61 

N 

2.10 

S 

1.81 
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Fig. 5. Variation of bed temperature (a), entropy generation (b) and second law 
efficiency (c) with different height ratios. 


{0) and two operational velocity values. In Fig. 5b, the accordance 
of the trend of the total entropy generation in lower operational 
velocity with the bed temperature trend of Fig. 5a is explained 
by the dominant effect of the heat transfer component of the entro¬ 
py generation. The effect of the fluid friction component of the en¬ 
tropy generation is negligible in lower operational velocities. As 
the operational velocity increases, the bed temperature reaches 
higher values due to decrease in bed-to-wall heat transfer coeffi- 


Fig.6. Combustor efficiency along the bed height at different heat transfer areas for 
<P = 0.3 (U 0 = 5.00 m/s). 


cient and less contact time in the bed [33,34]. This phenomenon 
also causes higher entropy generation values. As for higher opera¬ 
tional velocities, a definite increase in the total entropy generation 
is observed, which indicates that the fluid friction component of 
the entropy generation also effects together with the heat transfer 
component of the entropy generation. This situation is caused by 
the hydrodynamic behavior of CFBs. The combustor efficiency de¬ 
creases along the bed height as the operational velocity increases 
and its value at every level in the combustor decreases as the value 
of 0 decreases. For 0 = 0.6 value, the amount of decrease gets big¬ 
ger, regardless of the operational velocity value. Regardless of the 
positioning of the heat transfer surfaces in the combustor, the ef¬ 
fect of the fluid friction component of the entropy generation is 
more effective on the total entropy generation as the operational 
bed velocity increases. The difference observed in the value of 
the second law efficiency for different values of 0 indicates that 
the positioning of the heat transfer surfaces is of utmost impor¬ 
tance in the optimum performance of CFBs as shown in Fig. 5c. 

The effect of heat transfer tube area on combustor efficiency is 
also obtained at constant cooling water mass flow rate and at dif¬ 
ferent heat transfer area ratio, R h e which is the ratio of overall heat 
exchanger area to heat exchanger length in the bed (Fig. 6). During 
the simulations, bed is operated at fast fluidizing conditions and 
coal particle size is in the interval of 0.0125-6.5 mm in diameters. 
Combustor efficiency decrease along the bed height at constant 
heat transfer area ratio. The decrease in efficiency can be observed 
at every point along the bed height as the area ratio increases. This 
phenomenon is explained by the negative effect of the increase of 
R he on the bed hydrodynamics which is further justified by the 
dominant effect of the fluid friction component of the entropy gen¬ 
eration. As shown in Fig. 6, the increase of the second law effi¬ 
ciency with the area ratio indicates that the increase of the heat 
transfer surface area is one of the key parameters in the optimum 
performance design of CFBs. 

The effect of heat transfer surfaces on CFB combustor efficiency 
along the bed height at different volume ratios is shown in Fig. 7. 
The increase of the total entropy generation as the 7 increases indi¬ 
cates that the fluid friction and heat transfer components of the en¬ 
tropy generation are more effective on the total entropy generation 
as shown in Fig. 7b. The increase of the total entropy generation 
along the bed height for different values of the 7, indicates that 7 
affects the bed hydrodynamics negatively. In Fig. 7b, the obvious 
difference of the trend of the total entropy generation with the 
bed temperature trend of Fig. 7a is a further explanation of this 
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Dimensionless Bed Height, z/Hjj (-) 

Fig. 7. Variation of bed temperature (a), entropy generation (b), and second law 
efficiency (c) along the bed height with different heat exchanger volume ratios 
(Ho = 5.00 m/s). 


situation. Efficiency decreases with bed height at each y value, and 
its value at every level in the riser increases as the value of y in¬ 
creases up to y = 0.50. The decrement amount, especially at lower 
y values, is bigger due to the variation in temperature distribution 
within the combustor. As the heat transfer volume ratio increases, 
the amount of heat drawn from the bed also increases, resulting in 
lower bed temperatures and higher temperatures in the heat 



Volume Ratio, (-) 


Fig. 8. Effect of the volume ratio on second law efficiency at each CFB zone 
(Ho = 5.00 m/s). 


exchangers. At a higher volume ratio, the combustor efficiency va¬ 
lue stays constant due to the moderate combustor temperature 
which is caused by a smaller entropy generation resulting from 
the lower temperature distribution within the combustor. 

In CFB combustors, amount of heat transfer and efficiency con¬ 
siderably varies between the bottom and the upper zones due to 
different hydrodynamic characteristics of each zone. This variation 
is examined through the simulation results at the different volume 
ratio, for each zone (Fig. 8). At each zone as heat exchanger volume 
ratio increases, the average values of second law efficiency are get¬ 
ting higher values and it stays constant especially after y = 0.3 va¬ 
lue because of the lower temperature distribution within the 
combustor. 

Water mass flow rate and heat exchanger tube diameter are 
also important design parameters in CFB combustors. In Fig. 9, sim¬ 
ulation results are presented at two different cooling water mass 
flow rates and heat transfer tube diameters for a constant height 
ratio {0 = 0.3). Second law efficiency is considerably affected with 
the tube diameter due to higher heat transfer coefficient at the 
water side. Smaller heat transfer tube diameter also affects the voi- 
dage ratio and contact area of the bed material considerably in the 
CFB as shown in Fig. 9. In the smaller heat transfer tube diameters 
and the larger water mass flow rates, the amount of heat with¬ 
drawn increases and results in a positive effect on the second 
law efficiency. 

4. Conclusions 

The main contribution of this study to science domain and the 
designers lies in the fact that it proves that the heat exchanger sur¬ 
faces have an irrefutable importance on the quality of energy to be 
derived from fluidized bed combustion which also involves many 
complicated factors such as chemical reactions and physical trans¬ 
port processes. In the real applications, especially in real fluidized 
bed applications where the heat derived from the bed is of utmost 
importance; the design and positioning of the heat exchangers di¬ 
rectly affect the quality of the energy system. Thus the thermody¬ 
namic analysis of heat transfer surfaces in circulating fluidized 
beds is very important to show how the parameters affect the effi¬ 
ciency of a fluidized bed combustor. In this study, the thermody¬ 
namic second law analysis of heat transfer surfaces in CFBs is 
analyzed theoretically in order to define the parameters that affect 
the system efficiency using a previously developed 2D CFB model. 
Combustor efficiency is examined under various conditions and for 
different types of heat exchanger surfaces. As a result of this study, 
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Fig. 9. Variation of bed temperature (a), entropy generation (b) and second law 
efficiency (c) with different water flow rates (<P = 0.3). 

it may be claimed that heat transfer surfaces are extremely impor¬ 
tant on the system efficiency. It is observed that there is a general 
decrease tendency of the second law efficiency of the CFB combus¬ 
tor along the bed height and that the location of heat exchanger 
surfaces in the bottom zone increases the efficiency. As the opera¬ 
tional velocity decreases the efficiency of the system is also af¬ 
fected negatively. An increase of the R he -area ratio affects the 
system efficiency positively. When compared with the volume ra¬ 
tio (7) and area ratio (R he ), the effect of the height ratio ( 0 ) is ob¬ 
served to be dominant. Efficiency at every level in the combustor 
increases as the value of y increases up to 0.50. At higher volume 
ratios, combustor efficiency value stays constant. Water mass flow 
rate and heat exchanger tube diameter are also important design 
parameters in CFB combustors. Second law efficiency is consider¬ 


ably affected with the tube diameter due to higher heat transfer 
coefficient at the water side. Smaller tube diameter also affects 
the voidage ratio and contact area of the bed material considerably 
in the CFB. 
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